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ABSTRACT 

An open supersymmetric t-J chain with boundary fields is studied by means of the 
Bethe Ansatz. Ground state properties for the case of an almost half-filled band 
and a bulk magnetic field are determined. Boundary susceptibilities are calculated 
as functions of the boundary fields. The effects of the boundary on excitations 
are investigated by constructing the exact boundary S-matrix. From the analytic 
structure of the boundary S-matrices one deduces that boundary bound states are 
formed for sufficiently strong boundary fields. 
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I. INTRODUCTION 



Recently there has been renewed interest in one- dimensional impurity problems and the 
related problem of one-dimensional Luttinger liquids with boundaries |^,[|,f|,|j|||,[/]|| . The 
main focus of these investigations has been the effects of Kondo-like impurities and effects due 
to potential scattering in Luttinger liquids. These impurity problems are closely related to 
open 1-d systems with boundary fields. Some of these systems are integrable and can be solved 
exactly by Bethe Ansatz P, |IT||iT| , |T2"| , |il! , TA \ . In particular, in ]T[ an anisotropic Heisenberg chain 



with open boundary conditions was studied. It is the purpose of the present work extend the 
investigation of to the case of the t-J model, which is a Luttinger liquid with both spin 
and charge degrees of freedom. 



In |TH| a trigonometric generalization of the supersymmetric t-J model with open bound- 
aries was constructed by means of the Quantum- Inverse Scattering Method (see e.g. f25j). 
This generalized the previous work by Forster and Karowski on the quantum-group invariant 
case fll6 |. Here we study this model at the rational point, for which it reduces to the supersym- 
metric t-J model with open boundaries and boundary fields. The reason for this restriction 
is that the trigonometric model in general leads to a non-hermitian bulk hamiltonian. The 
only exception is the hyperbolic regime, but there the spin excitations are gapped and thus 
irrelevant for the low-energy physics of the model. The hamiltonian we consider in the grand 
canonical ensemble is given by 

H = -V (fj £ ct a c j+lja + c) +1>aCj ,)j V 

+2 J2 Sj ■ S j+1 - + J2 nj + n j+1 - HS* - fiN + H a(3 , (1) 

where V projects out double occupancies, Sj are spin operators at site j, rij = cj^c^j + cj jC J - ) ,[ p , 
and the four possible choices of boundary hamiltonians H a p compatible with integrability and 
conservation of total spin in z-direction and particle number are given by 

H aa = K n l + h L n L , H ab = h^S* - — ) + h' L Tl L 

H ba = h[ ni + h L (S z L - ^) , H bb = h^Sl - ^) + h L (S' L - ^) . (2) 

Here rij = 1 — rij^ — rij^ is the number operator for holes (unoccupied sites) at site j. To 
simplify the computations we constrain ourselves to the regions h e (0,2), h' 6 (0,1). It is 
straightforward to extend the analysis below to other ranges of the fields. For later convenience 
we define the quantities 

{2 — p- for aa, ba ( 2 — -M- for aa, ab 

, Sl=\ t ■ (3) 

1 — pj- for ab, bb [ 1 — ^- for ba, bb 

In what follows we always assume that S\ and Sl are noninteger numbers. We note that 
for zero boundary fields (0) exhibits a global s/(l|2) symmetry ||16|| . In the present work we 
perform a detailed study of the boundary effects in the model defined by ([j]), paying particular 
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attention to the influence of the boundary fields. After some technical preliminaries we turn to 
an analysis of the ground state properties. We find that the zero-temperature susceptibilities 
exhibit some interesting singularities, which we argue to be related to the formation of bound 
states near the boundaries. We then study the interaction of elementary excitations with 
the boundaries by computing the exact boundary S-matrices. We find that boundary bound 
states can be formed for sufficiently strong boundary fields. We concentrate on the case of 
band fillings close to one (corresponding to the particularly interesting case of the lightly doped 
Mott-Hubbard insulator) for which it is possible to obtain explicit analytical results. However 
it is straightforward to extend the present analysis to arbitrary band-fillings by solving the 
integral equations (|l4l ) numerically and then numerically integrating fljjl) . 

Taking the rational limit of the Bethe Ansatz equations derived in |L5| we obtain 



Nh+Ni, N h 

Wap(h) (ei(A fc )) 2L = Yl e 2(Afe-A i )e 2 (A A; + A j )ne_i(A fe -A / (1) )e_i(A fc + A / (1) ) 

3i=k 1=1 
N h +N L 

1 = C^(Af J ) n ex^-A.MA^+A,-), (4) 
i=i 

where e n (x) = — — ^ anc i a, P = a, b. The boundary terms are given by 

X 2 

?7 aa (A) = 1 , r] ab (\) = -e_ 5l (A) , i] ba (\) = -e_ Si (A) , r] bb (X) = r} ab (X)r} ba (X) 

C»(A) = 1, Cafc(A) = -e_ 5i (A) , C6a(A) = -e_ Sl (A) , ( aa (\) = ( a b(\)(ba{\) . (5) 

The restrictions imposed on h and h! are chosen such that in all these expressions the label x 
on e x (X) is positive with range (0, oo). The energy of a state corresponding to a solution of 
(HI) is (up to an overall constant, which we drop) 

N h +N x 1 H H 

E = E ij - £ T-^ + H(N l + N h ) + ^--)N h -(n+-)L , (6) 

3=1 4+ A J 1 Z 

where E aa = h[ + h' L , E bb = and so on. The reference state used to derive @ is the 

one with up-spin electrons at each site of the lattice. This leads to the constraint in (H) that 
the number iVj_ of down-spins must be smaller than or equal to the number of up-spins Np 
Solutions of (U) violating this constraint can lead to vanishing wave-functions and must be 
ignored. Eigenstates of (|1|) with iVj > iVj must be constructed by switching the reference state 
to the state with down-spin electrons at all sites. This leads to the same Bethe equations (f|) 
with Ni <-> iVj and different values for the quantities Sj 

{2 — y- for aa, ba ( 2 — S- for aa, ab 

, S L ={ L . (7) 

1 + y- for ab, bb [ 1 + y- for ba, bb 

Below we will mainly deal with situations for which iVj_ < Np However, when considering 
excitations over the antiferromagnetic ground states we will also consider the case > N^, 
for which we will make use of the procedure outlined above. 
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In order to simplify (f|) we make use of the 'string-hypothesis' 0, which states that for 
L — > oo all solutions are composed of real Aj^'s whereas the A's are distributed in the complex 
plane according to the description 

Kt = K + i(^--5) ,J = l-..n (8) 



where a = 1 . . . M n labels different 'strings' of length n. This string hypothesis is naturally 
identical to the one for the model with periodic boundary conditions. The imaginary parts of 
the A's can now be eliminated from (|J) via (H). Taking the logarithm of the resulting equations 
(for M n strings (§) of length n and N h A^'s (note that J2™=i nM n = N± + N h ) we arrive at 

2tv 1 \ n 1 

T 4 n = (2 + j)e{^) - - E OmniK - a™) + e mn (K + ^) 

(m/3) 

I N h yn _ ^(1) yn , ^(1) 

+7 E 9{p^=L-) + 6{^±^) + \$\K) , a = 1 . . . M„ 
Li -y n n ±j 

\j, = \Y. ^? - K) + o(\? + K) + 1^?) ,7=i- (9) 

(net) 

where 7™ and J 7 are integer numbers, 6{x) = 2 arctan(2a;), 

O n ,m(x) = (1 - W^Cl—^— |) + 2 ^ + • • • + 2 0( - - - ) + , 

|n — m| |n — m| + 2 n + m — 2 n + m 

(10) 

and the boundary contributions are given by 



K 



in) 
ab 



(A)=E^( Zi 9/ <? )» ^( A ) = ^(^) 

j— [ n + l-2/-6i -6 L 

A , A x 



4? (A) =E^ L, 1 „, c ) » = e ( 

i=i 



^(A) = ^ J (A) + ^(A), ^(A) = 0, 
«W(A) = 0, ^ aa (A) = a; 6a (A) + ^ ab (A) . (11) 

The ranges of integers J™ and J 7 are 

oo 

/" = 1,2,... L + M n -2 E mm{m,n}M m + N h , J 7 = 1, 2, . . .i^ + JV h - 1 . (12) 

m=l 

There are two differences as compared to the case of periodic boundary conditions 



I7| . p^ , p!9l , p0| , ^Tl , |22| , ^3| , p4| : first there are additional terms (|TT|) , and secondly the integers I™ 



^ As far as the present work is concerned we do not need to explicitly consider complex solutions of 
the Bethe equations and all of our results are really independent of the string hypothesis. 
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and J 7 take different values. The allowed range of the integers J™ and J 7 reflects the fact that all 
solutions of (|4]) with one or more roots Xj or A^ having vanishing real parts must be excluded 
as they lead to vanishing wave-functions. This restriction leads to constraints on the allowed 
values of the integers J™ and J 7 : the 7™ range from 1 to L + M n — 2 J2m=i mm { m , n}M m + N h , 
the solution with J™ = being excluded. Similarly J 7 range from 1 to N± + Nh — 1 and is 
again excluded. 

For zero boundary fields (k\j = 0, u>y = 0) we can costruct a complete set of 3 L states from 
the Bethe Ansatz states defined in the above way: the model (|l|) with vanishing boundary fields 
is sZ(l|2)-invariant and all Bethe states are highest weight states of s/(l|2) P^ITSI ]. Additional 
linearly independent eigenstates of (HD can be constructed by acting with the sZ(l|2) lowering 
operators on the highest- weight states. The total number of states obtained in this way is 3 L as 
can be proved in the same way as for the periodic t-J chain [|17]] (the necessary combinatorics 
are identical). Thus we obtain a complete set of eigenstates of ([[]). 

For nonvanishing boundary fields the situation is more complicated as the s/(l|2) symmetry 
is broken by the boundary conditions. Therefore we cannot use the symmetry generators 
to construct additional states from the Bethe Ansatz states and are left with the a priori 
incomplete set of eigenstates given by (|9|) and ([12]). For the present purposes this is inessential: 
the ground state is always a Bethe Ansatz state, as are the states needed to extract the 
boundary S-matrices. Ground state and excitations can be constructed from (||) in a standard 
way (see e.g. |2~5||). The ground state is obtained by filling all allowed vacancies of integers 
and J 7 up to maximal values J max and J max , which corresponds to filling two Fermi seas of 
rapidities A„ between and A and A^ between and B. The actual values of A, B (and thus 
J max and J max ) depend on H and \i and are determined below. We are interested in the case 
of a small magnetic field H and a close to half-filled band (/i ~ 2 ln(2)), for which A ^> 1 and 
Bel. As is shown in Appendix A the ground state energy per site (for the four possible 
sets of boundary fields) below half-filling is given by 

E -^ l ~ HSZ = ec (0) - 2a* + -L £ dX e s (A)4(A) + -L £ dX e c (A)o/ y (A) 

-^M0) + M - f - 2E ij ] + oiL- 1 ) , (13) 

where i,j — a,b and where the dressed energies £ C (A) and £ S (A) i are given in terms of the 
coupled integral equations 

/•A rB 

£ S (X) = — 27rai(A) + H — / dfi a 2 (X - y) eJ/j,) + dfi ai(A — n) e c {n) , 

J-A J-B 

e c (X)=fi- — + dy a x (A - y) e 8 {ji) . (14) 

Z J-A 

Here a n (X) = tt- n n $ . For later use we define 



i These can be shown to be (minus) the energies of the order one contributions to the elementary 
charge and spin excitations |22]. 
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G X (X) 



2tt 



du exp(— iu)\) 



exp(— x\^\) 1 

2cOsh(|) ~ 7T 



Re [P( 



1 + x 



+ *A) , 



(15) 



where x is real and where (3{z) = \ 
asymptotic behaviour of G X (X) for 



— t/>(~ ) . Here ip(z) is the digamma function. The 
arge / ^> 1 and A 3> x is 



Below we will also need the small- A asymptotics of G?i(A), which is given by 



(16) 



d(A) = ^[21n(2) + 2 E(-l) n (l " 2- 2 ")C(2n + 1)A 

Z7r n=l 



2in 



|A| < 1 



(17) 



II. WIENER-HOPF ANALYSIS FOR THE DRESSED ENERGIES 
In this section we analyze the coupled integral equations (|T4|) by means of Wiener-Hopf 



techniques |30| (for detailed expositions see e.g. f31| , |2"9|j ). As ( |i4|) are similar to the analogous 
equations for the densities in the periodic t-J chain the necessary steps are the same as in 
|19|| . However as we will need more explicit answers than are given in |19| for determining the 



boundary contribution to the ground state energy we briefly summarize the most important 
steps below. After Fourier-transforming, the first equation of (LI) can be turned into a Wiener- 
Hopf equation for y(X) = e s (X + A) 

y(X) = -2nG (X + A) + ^+ f° du [G x {\ -v) + G X (X + v + 2A)\ y{v) + CG (X + A) , 

z Jo 



where C = J^ B du exp(7ri/) e c {y). Here we have used the fact that A ^> 1 and B < 1 to 
approximate §^ B dv G (X — v + A) e c {y) ~ Go(A + A) J^ B du exp(7rz/) e c {y). The quantity 
C is determined self-consistently below. Eqn fll8|) can now be solved by iteration y(X) = 
?/i (A) +y 2 (X) + where 

yi (X) = -2nG (X + A) + ^ + [°° du G X (X - v) y x {v) + CG (X + A) , (19) 

2 Jo 

poo poo 

y 2 (X) = dvG 1 (X + v + 2A) y x (u) + dvG x {X- v) y 2 {v) . (20) 
Jo Jo 

These equations can be solved in a standard way through a Wiener-Hopf factorization. The 
result for y x is obtained in complete analogy with e.g. the Appendix of [^|] 

yt(u) = G» l iHG ~(°) * - ^(2n - C)G~ (-in f^^ } + 0(exp(-2vrA)) . (21) 

Here the Fourier transform yf{uj) = Jq°° dX y(X) exp(iXu>) is analytic in the upper half-plane 
and G ± (lv) are analytic functions in the upper/lower half-plane factorizing the kernel 1 + 
exp(-M) = G+(u)G-(u) 
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The equation for ^(A) is more difficult to solve. They key is to use the fact that X+X'+2A ^> 1. 
Using the asymptotic behaviour (ITH) of Gi(A) in the expression for the driving term D(X) = 
Jq 00 dX' G\(X + A' + 2A)yi(X') and then performing a Laplace transformation we obtain 

^(A) - — / dx exp(-2Ax)exp(-|A|a;)y 1 + (zx)[x + — + •••] , (23) 

4:71 JO 12 

where the expansion in x corresponds to the asymptotic expansion of Gi(A + A' + 2A). It 
is clear that due to the strongly decaying factor exp(— 2Ax) the leading contribution to the 
integral comes from the small-x region. Inserting the expression fl2"5| ) for the driving term into 
and then following through the same steps as in the analysis for y\ we arrive at 

* f°° J™ o A„\\„ , x3 _| j^Nsii"^ 1 ) 



yjfo;) ~ G + (c<j)— / dx exp(— 2Ax)\x + 

4:71 Jo 



12 to + ix 



„ G+H^H fdx exp(-2^) 1 + ^' nW + ---. (24) 

By means of a similar analysis further corrections to y{X) can be determined. As far as the 
physical quantities determined below are concerned y 3 , y± etc give rise to contributions much 
smaller than those due to yi and y%. We are now in a position to determine the limit of 
integration A as a function of the magnetic field. By definition e a (±A) = = y(0), which 
leads to 



A = j5 W + N£^-£)) i + 

71 71 A7Tln{H) 

Using fl2l|) and ( f24"D we can now solve the integral equation ( |i~4"l) for e c (A) 

£ C (A) = e (A) + f B duldiX -u) + 2QCQ Sh(7rA) exp(rru)}e c (u) , (26) 

J-B 271 

where a = - exp(— 2ttA), £o(A) = \i — 2txG\{X) — 2acosh(7rA), and where we have neglected 
terms of order o(exp(— 2txA)). Here the term proportional to a originates in the C-term in (pl|). 
Equation fl26|) can now be solved by iteration as 73 1 (corresponding to \p\ = \2 ln(2) — fj,\ <C 
1) with the result 

e c (A) = fj, - (2tt + g)G 1 (X) - 2a cosh(TrA) + Ofjm) + 0(/2 2 ) + 0(a 2 ) , (27) 

where 

/ ^' 1 -777^7K + ---) > ^ = |-r==02 + 2o)*. (28) 



8tt^ 21n(#) 3 y ^) 
The boundary of integration B defined via e c (±B) = in this order is given by 



3C(3) 



81n(2) 1 , 

// + 2a + =(/Z + 2a )'■ 

3vr ./6C(3) 



(29) 



C is determined self-consistently to be — 2((3)B 3 . The higher order (in B) contributions to e c 
and B do not contribute to the singularities in the thermodynamic quantities and therefore 
have been dropped. 
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III. GROUND STATE PROPERTIES 



We are now in a position to determine bulk and boundary contributions to the energy (|13|). 
The bulk energy per site is found to be 

e bu ik = - [p. + 2 ln(2) + 2a + 2 ln(2) ^B 3 } , (30) 

7T 

from which we can determine the leading contributions to the zero-temperature magnetization 
per site, magnetic susceptibility, density and compressibility close to half-filling in a weak 
magnetic field 

debuik H . I 
mbuik — tttt — tt^U ~~ 



dH 2tt 2V 2\n.{H)' 

1 1 21n(2) H 2 

XH )b u lk - — 2 (l ~ ^y) + -^4^(1 - ^H)\/ 6am + 2a) 4 




n de bu ik _ 2 ln(2) 

Ubulk — — 1 



2(/2 + 2a) ^ 



d/M vr \\ 3C(3) 

1 aD^fc 21n(2) 1 
D buik 9fM n jQ((3)[fi + 2a} 



in agreement with the expressions for periodic boundary conditions P^| , pO| . We note that both 
the magnetic susceptibility and the compressibility diverge when we approach half-filling. 



Contributions to the surface energy, i.e. all terms proportional to L^ 1 in (13), can be 
divided into boundary- field dependent ones E^ a ^ and contributions due to the "geometry" 
i.e. openess of the chain E° so that we can write for the four permitted sets of boundary 
conditions 

E boundary = E° + , a,/3 = a,b . (32) 
The boundary field independent contributions are easily determined 



--4{-^( 1+ w) + '*--& +2o)l+ -]- < 33 ' 

We note that for zero bulk magnetic field H = and half-filling fi = 2 In 2 we obtain E° = 
f — In 2, which is the correct result for the surface energy of the open XXX Heisenberg chain 



I3lj35| . By differentiating the surface energy with respect to the thermodynamic parameters 
H and fi we can evaluate the surface contributions to particle number and magnetization in 
analogy with e.g. the treatment of the Kondo model |31| (see also |L|). It is reasonable to 
assume that these contributions are concentrated in the boundary regions, e.g. we interpret 
the surface contribution to the particle number to lead to a depletion/increase of electrons in 
the "vicinity" of the boundaries. 

The leading contributions to the boundary magnetization, particle number and suscepti- 
bilities due to E° are 
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M o = 1 L ln(ln(ff)) \ o _ 1 / ln(ln(g)) \ 

' "'" 41n(#) V 21n(#) )^'~> AH(\n(H)) 2 \ 2\n(H) 

n °=^ i+ \\^¥)' *°=- / 4I " 2+T • < 34 » 

\ 3 C(3) / 27r A /6C(3)(^ + 2a) 



We first note that boundary region exhibits a stronger magnetization as compared to the bulk, 
i.e. M ° = —TrwT—n, which is much larger than one for the small magnetic fields considered 

m bulk 2H\nH^ ° ° 

here. The boundary magnetic susceptibility is seen to diverge for H — > 0. Following ]I[ we 
interpret this as an indication for the presence of a magnetic bound state in the boundary 
region. The magnetic behaviour is similar to the one for the XXZ spin chain with an open 
boundary studied in 

The leading contribution to the boundary particle number is |. Because of the constraint 
of at most single occupancy at any site this increase in particle number (recall that we are very 
close to half-filling) must be spread out over large regions neighbouring the boundaries. This 
indicates that boundary effects spread deeply into the bulk. The boundary compressibility 
for zero boundary fields is seen to be negative and to diverge as we approach half-filling. The 
type of singularity is the same as for the bulk. Combining the results for magnetization and 
particle number we see that there is a tendency for spin-up electrons to get pushed towards 
the boundary. 

The leading order boundary- field dependent contributions are expressed (see fll3|)) 

in terms of the quantities 

rA rB 

e b (S) = dv a s (y) e s {v) , e a (S) = dv a s {y) e c (v) . (35) 

J-A J-B 

where according to (|T3| ) for the four possible sets of boundary conditions 

£ (Q/3) = \ M-St) + e a (-S L )] + E a(3 , a,/3 = a,b. (36) 

The leading contribution to the quantity e a (S) can be easily determined for the case S ^> 
B, in which we can expand as{v) in a power series in v and then perform the elementary 
integrations using the expression (B?j) for e c (y). For S B we instead expand e c (v) in an 
infinite power series (using that G\{y) is a smooth function around zero), then perform the 
integrations, resum the result, and then retain only the leading terms. This results in 



i + if g >> B 

rjS)={ * S (37) 
£ c (0) + ||(/2 + 2a + ...) + -.. if S<.B. 

The analogous computations for €b(S) are more involved, so that we give a brief summary of 
the necessary steps in Appendix B. We find the following result for the leading behaviour 

e b (S) = -G 5 (0)(27r + 2C(3)5 3 ) + -+ HH) . (38) 

2 {-T-^s H HH) + ... ifS»A 
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A. Contributions due to small boundary fields of type a 



The contribution to the boundary energy is given by |e a (| — 2) with | — 2 ^> B = J 2 ^||^ ■ 
(This defines what we mean with "small" boundary field). Here h is a boundary chemical 
potential. We define the quantity a = TZ^^f^^ which in the present case is much smaller 
than one. We obtain the following contributions to boundary magnetization/particle number 
and susceptibilities 

X " e = -(p + 2a)- 1 . (39) 

7T 

We see that a small boundary chemical potential leads essentially to the same type of diver- 
gences as are present in the bulk. As expected electrons get pushed away from the boundary 
although the effect is small. By differentiation with respect to the boundary chemical potential 
we can evaluate the average number of electrons at the boundary site 

<*»> - 1 - 2C f (1 B _ 3 hy . ( 4 °) 

where B is given by (^9|). We see that the electron number is larger than the bulk value. 
This is consistent with the above observation that an open boundary without field leads to an 
increase in the elctron density in the boundary region. 



B. Contributions due to large boundary fields of type a 
Here the boundary chemical potential is taken large, by which we mean that < f-2 < B. 



We again use the notation a = j^JqtM, but now a ^> 1. We find 



2(/2 + 2a)^ 1 



6C(3) : 

31n(2) 1 , , , 

7rJ6C(3)(/i + 2a) 4 ™ 



The magnetization is again proportional to H and the magnetic susceptibility can only diverge 
at half-filling. The large boundary field yields however a contribution of — | to the boundary 
particle number, which indicates a strong depletion of electrons in the boundary region. This 
is in accordance with the expectation based on a naive analysis of the hamiltonian (P that 
large boundary chemical potentials (with our choice of sign in (Jl])) should favour the presence 



9 



of holes in the boundary region. The compressibility exhibits a stronger divergence than the 
bulk if we approach half-filling keeping a fixed. The average electron number at the boundary 
site is found to be 

(n e ) = l-^-jQam + 2a) , /wl, (42) 
which is less than the bulk value. 



C. Contributions due to large boundary fields of type b 



Let us first consider the region where h\ In if | ^> 2tt, which corresponds to the case 5 C i. 
By straightforward differentiation we find 



1 



M° = — 



h-l 



HB 



4 2/i In if 4tt 



2h ' 



^ K 2W 



1-h 1 

+ 



N = — 



2hH\n 2 H ' 4tt 3 
1 



VK 2h ' VK 2h' 



B 



ff 2 



6vr 2 C(3) J B / 



1 

1 

2~k 



2h ' ^ K 2h> 



'^(— )-^(-) 
VK 2h ' ^ K 2h } 



2(fl + 2a) 
\ 3C(3) 
1 



^6C(3)(/i + 2a) ' 



(43) 



Thus the boundary field contributes to the singularity in the magnetic susceptibility for large 
boundary fields h. Furthermore there is a constant contribution — | to the boundary magneti- 
zation, which indicates a surplus of spin-down electrons in the boundary region. The negative 
sign stems from the fact that the boundary field is antiparallel to the bulk magnetic field. 
The boundary particle number contribution is always small (and leads to a depletion of the 
electron number in the boundary region) and the compressibility exhibits the same type of 
divergence as the bulk. By differentiating with respect to the boundary field we can calculate 
the (S z — ?y) at site 1 (L) for ab (ba) boundary conditions. The result is 



(S z - 



n 



1 r c 

2 ~ Jo 



duj uj 



H 



exp(— u) 
1 + exp(-hu) ~ 2h 2 In if 



(44) 



To get a rough idea of the contribution 
1 and h = 2, which are fj and 0.91597 



This is always finite in the range of h considered, 
of the integral we note its respective values for h ■ 

(Catalan's constant). The contribution proportional to H is always small. The susceptibility 
is finite, which means that the spins and charges at the boundary itself do not contribute to 
the bound state responsible for the singularities in the susceptibilities. 



D. Contributions due to small boundary fields of type b 
Let us now turn to the region h\ In if | <^ 2n of vanishingly small boundary fields. We find 
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M b = , x°h 



hln(H) b h 



2n 2 H 



3C(3) 27r v /6C(3)(/z + 2a) 



Note that we cannot take H — > without taking the boundary field to zero first. Thus the 
magnetization is always small. However the boundary magnetic susceptibility may or may 
not diverge for H — > 0, depending on how fast we take take the boundary field to zero as 
compared to the bulk field. The result for (S z — ! |-) is now found to be 

(s . _ £ = _™ ^ (46) 

This is again small and vanishes for half-filling and zero bulk field in accordance with ffl. The 
corresponding susceptibility is again finite and the spins/charges at the boundary site do not 
contribute to the bound state. 



IV. BOUNDARY S-MATRIX 

In this section we study the effects of the boundary on gapless excitations. For simplicity 
we set the bulk magnetic field H to zero. As the elementary excitations in the bulk are the 
same for the periodic and the open chain we begin by giving a thorough discussion of the 
interpretation of the spectrum in terms of elementary excitations for the periodic system. 
After reviewing the known results of we present a conjecture concerning the sZ(l|2) 

descendants of the holon and spinon states obtained from the Bethe Ansatz. 

There are two kinds of gapless elementary excitations in the supersymmetric t-J model, 
associated with spin- and charge degrees of freedom respectively. For the periodic chain they 



have been extensively studied in fl2l|j22|l (see below for a discussion of the sl(l\2) structure of the 
excitations). The spin excitations are called spinons and carry spin ±| and zero charge. They 



are very similar to the spin- waves in the Heisenberg XXX chain [38| . The charge excitations are 
called holons and antiholons, carry zero spin and charge =Fe. Holons correspond to "particles" 
in the charge Fermi sea of A^'s and are thus associated with a physical hole, whereas antiholons 
correspond to "holes" (unoccupied A^'s) in the charge Fermi sea. At half-filling only holons 
can be excited as the charge Fermi sea is completely empty. The excitation energies are given 



by £ S)C (A) (|14D , whereas their physical momenta (for the periodic chain) are given in terms of 



the quantities p S)C (A), which are solutions of the integral equations 

pA i-B 

p s (A) = -0(A) - / dv a 2 (X - f)pa(f) + / dv aA\ - v) p c (i/) 

J -A J-B 



A 



Pc (A) = / dv ai(A - v) p». (47) 



A 



The mometum of e.g. a holon-antiholon excitation is given by P cS = p c (A p ) — p c (A h ) where A p 
and A h are the spectral parameters of the holon and antiholon respectively. We thus would 
define the physical holon momentum as p c (A p ) = p c (A p ) — p c (B). At half-filling the spinon 
(p s ) and holon (p c ) momenta are given by 
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Ps(X) = 2 arctan(exp(7rA)) — n , 

Pc(A) - 2 + « ln 1^(^)1X1- f) J • (48) 

The "order one" contributions to the spinon and holon energies at half-filling take the simple 
form e s (A) = 27rG (A) and £ C (A) = 21n(2) — 2nGi(\). The spinon dispersion is thus of the 
Faddeev-Takhtajan form £ s (p) = 7rsin(p). The holon dispersion cannot be written in closed 
form so easily. 

So far we have not discussed the role played by the s/(l|2) symmetry in the classification 



of eigenstates. As was shown in |17j all eigenstates of the hamiltonian obtained from the 
Bethe Ansatz are highest weight states of the s/(l|2) symmetry algebra of the model. Addi- 
tional eigenstates of the hamiltonian can be obtained by acting with the s/(l|2) generators 
(recall that we are still discussing the periodic chain). This leads to the picture put forward 
in [|T7| for the structure of excitations over the antiferromagnetic ground state, i.e. spinon 
and holon/antiholon excitations are really only the highest weight states of s/(l|2) multiplets. 
However all the additional excitations constructed by acting with the symmetry generators 
with the exception of the spin-S77(2) descendants can be easily shown to have a gap pro- 



portional to the chemical potential (see also |26|]). Therefore we can obtain a complete set 



of gapless excitations by taking into account spinon and holon states plus their spin-SU(2) 
descendants. However, one can furthermore argue that in the thermodynamic limit (i.e. if we 
neglect all finite-size corrections) actually all the gapped s/(l|2) descendants obtained by act- 
ing with the symmetry generators can be viewed as being incorporated in multiparticle spinon 
and holon/antiholon excitations. Let us first discuss the situation for a large finite chain. Here 
a complete set of states is given by first finding all sets of spectral parameters solving the Bethe 
equations (f|). Each such solution yields the wave-function of a highest- weight state of s/(l|2) 
with a given fixed momentum, and a complete set of states is obtained by taking into account 
the s/(l|2) descendants of the highest-weight state. As we approach the thermodynamic limit 
we identify one-parameter families of highest weight states as elementary excitations, the free 
parameter being the rapidity (which is directly related to the momentum) of the particle. Thus 
in the thermodynamic limit we are interested in multiparameter families of excited states and 
the strict counting of states possible in the finite volume loses its meaning. This is the reason 
why the sZ(l|2) descendants of the spinon and holon/antiholon excitations do not have to 
be taken into account separately anymore in the thermodynamic limit. From the Algebraic 
Bethe Ansatz construction we know that the symmetry generators can be represented as the 
infinite spectral parameter limits of the elements of the monodromy matrix [I7]f27]]. On the 
level of the Bethe Ansatz states this means that the action of the symmetry generators can be 
implemented by taking a spectral parameter in an appropriate Bethe Ansatz state to infinity. 
If we therefore take k rapidities of an n parameter family of excited (highest weight) states to 
infinity we produce an n — k parameter family of exited states made of s/(l|2) descendants! 
This means that the family of s/(l|2) descendant states does not have to be taken into account 
separately anymore, as can equally well be interpreted as "sitting on the boundary" of the 
n-parametric family of highest-weight states. In this way we obtain an equivalent yet different 



"quasiparticle-interpretation" in the spirit of McCoy et. al. |33 . 

As the simplest example let us consider the s/(l|2) descendants of the antiferromagnetic 
ground state at half-filling, which sits in a multiplet of dimension four, the other three states 
being Q a \GS) (a =t, I) with momentum zero and energy E = 21n2 and QiQ^\GS) with 
momentum zero and E = 4 In 2. The state QjjGS) can be obtained from the spinon-holon 
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scattering state (or more precisely the two-parametric family of states) by taking the spectral 
parameters of both the holon (A p ) and the spinon (X h ) to infinity: indeed the quantum numbers 
match and e s (X h ) + e c (A p ) -> 21n2 and p s (\ h ) + p c (A p ) -> 0. The state Q T |GS) can then 
be obtained by acting with the spin raising operator . Similarly the state Qj.Qt|GS) is 
obtained from the (two-parametric) holon-holon scattering state by again taking both spectral 
parameters to infinity. 

On the basis of the above discussion we therefore propose that the quasiparticle interpre- 
tation of the order one excitation spectrum in terms of two spinons with spin ±| and holon 
and antiholon excitations put forward in [E2J does actually already incorporate the complete 



s/(l|2) structure. For the half-filled band it is straightforward to show El by using the distribu- 
tion of integers ( |12|) and the highest weight theorem that all gapless excitations are scattering 
states of two spinons and one holon with the superselection rule that the number of spinons 
plus the number of holons is even. Thus the excitation spectrum of the half-filled t-J model 
is described by a 577(2) x U(l) scattering theory. 



The scattering matrix has been determined by means of Korepin's method f34]j25| in pl] , p2| . 
At half-filling the spinon-spinon S-matrix S(X) and the spinon-holon (sc) and holon-holon (cc) 
scattering phases are given by 




ie 7TX ' 



(49) 



where I and P are the 4x4 identity and permutation matrices respectively. Below half-filling 
the S-matrices are given in terms of the solution of integral equations. This concludes our 
discussion of the excitation spectrum of the periodic t-J model. For the case of the open chain 
the elementary excitations are identical as are the bulk S-matrices. What remains in order 
to completely specify the scattering of elementary excitations is to determine the phase-shifts 
acquired by the spinons and holons when reflecting from one of the boundaries. Note that due 
to the particular form of the boundary interactions (no particle production or transmutation) 
it is clear that the boundary S-matrices for holons and spinons are diagonal and thus reduce 
to phase-shifts for the physical states. 



A. Boundary S-matrix for the exactly Half-Filled Band 

For the case of the exactly half-filled band the boundary S-matrix can be determined by 
the method introduced by the Miami group in for the case of the spin | XXX Heisenberg 
chain, which generalizes the methods of Korepin |53] and Andrei et al [3"6,3"7|]. An alternative 



method was introduced in yj] and can be seen to lead to the same results. In order to determine 
the boundary phase shifts for spinons and holons we will study the spinon-holon scattering 
state. The method of |35[ is based on the following quantization condition for factorized 
scattering of two particles with rapidities Ai_2 on a line of length L (see also [[|) 



§The computation is analogous to the appendix of [32 
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exp(2iLp(Ai))5(Ai - X 2 )K 1 {X 1 , h 1 )S{X 1 + X 2 )K 1 {X 1 , h L ) = 1 , (50) 



where p(X) is the expression for the physical momentum of the corresponding (infinite) pe- 
riodic system, 5(A) are the bulk scattering matrices for scattering of particles 1 and 2, and 
Ki(X, h) are the .fT-matrices describing scattering of particle 1 with rapidity A off a boundary 
with boundary field h. For the present case of spinon-holon scattering with boundary fields 
preserving the two U(l) symmetries of spinon and holon numbers this condition turns into 
scalar equations for the scattering phases, which after taking the logarithm read 



2Lp s (X h ) + (p sc (X h - A p ) + ^j s (X\ ht) + <p(X h + A p ) + ij> a (\ h , h L ) = mod 2tt . 
2L Pc (A p ) + Vsc (A p - X h ) + ^ c (A p , h) + <p(A p + X h ) + ^ C (A P , h L ) = mod 2tt 



(51) 



Here X h and A p are the rapidities of the spinon and holon respectively. Comparing these 
conditions with certain quantities ( "counting functions" ) that can be calculated from the Bethe 

m 



Let us now discuss 

L 



Ansatz solution one can then read off the boundary phase-shifts ift s> 

this program for the spinon-holon scattering state characterized by choosing Mi = ^, = 1 
in the Bethe equations (|9|) . There are ■§ + 1 vacancies for the integers l\ and thus one hole in 
the Fermi sea of A*. We denote the rapidity corresponding to this hole by X h . The rapidity 
corresponding to the holon is denoted by A p . The Bethe equations read 



r+l 



§4 = (2 + \)6{x a ) - \ E e^i) + e(^±^] 



+j^ K if(X a ) + 



i 

L 



■ X a — X h X a + x h , 

. n > + U \ o 



9{X a - A p ) + 6{X a + A p ) 



2tT 1 2 I 1 r 

-t j =tY1 e ( AP - A «) + 6 ( AP + A -) + t^( ap ) - t K ap - ^ ) + ( ap + 

a=l 



(52) 



We note that the ground state at half-filling is identical to the one of the XXX Heisenberg 
chain and is obtained by filling the rapidities A^ between — oo and oo. In the limit L — > oo 
the distribution of roots X a is described by a single integral equation for the density of roots 
p s (X), which is of the same structure as (3.28) of fL3f . The main complication is that we need 



to take into account all contribution to order ^ and thus must deal with the fact that the roots 
are distributed not between — oo and oo but between two finite, L-dependent values —A and 
A. The integral equation is of Wiener-Hopf form but cannot be solved in a form sufficiently 
explicit for the purpose of determining the boundary phase-shifts. Following and we 
make the assumption that the contributions due to the shift of integration boundaries will be 
of higher order in i as far as the boundary phase-shifts are concerned and take A = oo (the 
validity of this assumption is discussed at the end of the section). The integral equation then 
can be solved by Fourier transform with the result 



Ps(u) = 2G (u) + - {Gi(w)[l + 2cos(cuA h )] + G (u)[l + 2cos(cuA p )] + / fi (w)} 



(53) 



where G x (u) 



exp(-f H) 
2cosh(f ) 



and where the contributions due to the boundary fields are 



fab(u) = G-ts^u) , fba(v) = G_i_ Si (^) , f bb (uj) = f ab (uj) + f ba (uj) , f aa (uj) = . (54) 
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For the further analysis it is convenient to define counting functions z s (\) and z c (X) 
z,(\) 



Z7T Z7T p = i 2 2 



,\ + \ h 



+ 0{\ - A p ) + 9(X + A p ) 



-+i 



z c {A) = — £ #(A - A Q ) + 9(A + A Q ) + -U y (A) - [fl(A - A h ) + 9(A + A" 

Z7T — , 111 Z7T L 

a=l 



(55) 



Note that for any root e.g. X a of (|52l) the counting function takes the integer value z s (X a ) = 
I a by construction. In the thermodynamic limit 4 times the derivative of £ S (A) yields the 
distribution function of rapidities p s (A). Straightforward integration of the density p s (A) 
yields the following results for the counting functions in the thermodynamic limit evaluated 
at the rapidities of the spinon and holon respectively 



2nz s {X h ) = 2L Ps (X h ) + ifsc{X h - A p ) + <p sc (X h + A p ) + s (A h ) = mod 2vr , 
-2nz c (A p ) = 2Lp c (A p ) + ip sc (A p - X h ) + <p sc {X h + A p ) + C (A P ) = mod 2tt , 



(56) 



where p s ,c{X) are the spinon/holon momenta ([48]), (p S c{X) is the bulk phase-shift for spinon- 
holon scattering (|49|), and 

roc (J^q _ _ „ , 

4> S {X) = - I — Gx(u)(l + exp(-iuX)) + G (u) + fij(w) exp(-iuX) , 



0c(A) 



00 du 



oo ILO . 



Gi(o;)(l + exp(-z^A)) - G {u) + fij(u) exp(-|-|) 



exp(— iuX) — uJij(X) 



(57) 



The last equalities in (|5q) hold due to the fact that evaluation of the counting function at a 
root of the Bethe equations yields an integer number. From these equations we can now infer 
the boundary phase shifts by comparing them with the quantization condition fl5ip . 

The scattering of spinons on a 6-type boundary with boundary field h is identical to the 
one in the XXX spin chain and the results are the same as |35] : the phase for a spinon with 
spin-up and rapidity A is given by 



r( 



2 J L \ L -r 2 ) L \ 4h 2 I 1 V Ah ^ 2 



r(i + f)r(i 



iA y( 2— h _i_ iX\ y( 2+h iX 

2 ' \ 4h 2' \ 4/i 2 



(5? 



The analogous phase for a spin-down spinon can be obtained in the following way [3l|: As 
pointed out above switching to the reference state with all spins down leads to a redefinition 
of the quantities Sj. The excitation constructed in a way analogous to the one above over this 
reference state has spin quantum number S z = — ~. In order to study negative boundary fields 
we need to keep track of the modification in the quantities «Si,l, which are now always positive. 
Repeating the above analysis for this case we obtain the following boundary phase-shifts for 
spin-down spinons 



X + i 



2-h 
2h 



A 1 2h 



(59) 
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The phases for scattering of spinons on a-type boundaries are very similar, e.g. 



r(|-f)r(i + f) 
r(i + f)r(i-f) ' 



(60) 



These expressions can be obtained from the o-type phases by setting the boundary fields 
to zero. Physically this means that the spinons do not "see" the a-type boundary fields at 
half-filling. 

The phases and ip^ accumulated by a holon scattering off a boundary of type a or b 
are given by 



Thus scattering of holons off the boundaries is influenced by both types of boundary fields. Let 
us now investigate the analytic structure of the above phase-shifts. In the physical strip for 
the spinon rapidity was defined by the condition |Im(A)| < 1 on the basis of periodicity of the 
expressions for momentum and energy. We propose the further constraint on the physical sheet 
that the imaginary part of the spinon momentum has to be positive (interpreting the spinons 
as particles). This implies that the spectral parameters should lie in the strip < Im(A) < 1. 
It now can be seen that the spinon boundary S'-matrices have no poles on this strip (note that 
for e*^ si there is no pole at A = i^r)- Therefore there are no boundary bound states in the 
region of boundary fields we consider here. Let us now turn to the holon boundary S'-matrices. 
The physical strip for the holon rapidity (for vanishing real part) is given by —1 < Im(A) < 0. 
We see that the holon boundary phase-shifts do not have poles on the physical sheet for the 
range of boundary fields considered here {h G (0,2) and h' G (0,1)). However we see that 
the pole A = ijZh °f dUD f° r a ~type boundaries starts crossing over onto the physical sheet 
for hi — > 1. This indicates that for h' > 1 a holon bound state forms. Indeed we find that 
for h > 1 a solution of @ exists where one of the roots A^ takes the value A = ijz^r- For 
h' > 1 this root is present in the ground state so that it indeed corresponds to a boundary 
bound state. We see that a sufficiently strong boundary field is needed for a bound state to 
be formed. 

Let us now dicuss the validity of our "incomplete" \ analysis of the densities/counting 
functions. As pointed out above a complete analytical treatment encounters significant tech- 
nical difficulties. However our results can be checked numerically in the following way: if the 
shift in integration boundary would lead to an additional term of order r, (|56|) should be 
incorrect to order one. We performed an explicit numerical evaluation of A p and X h corre- 
sponding to integers L/n and L/m with n and m fixed by solving the Bethe equations fl52|) 
for chains up to 600 sites. Through insertion of the numerical values of X h and A p into fl56|) 
it is then possible to check whether the neglected finite-size effects contribute to order one 
in the counting functions. We found no evidence for any missing contribution to the bound- 
ary phase-shifts. We therefore conclude that for the half-filled band the shift in integration 
boundary can indeed be neglected. 




(61) 
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B. Boundary S-matrix for the less than Half-Filled Band 

Let us now turn to the case of the less than half-filled band. It is a well-known fact that 
below half-filling the scattering matrices are not functions of the differences of the rapidities 
of the scattering particles any longer. Therefore we need to replace (|50|) by an appropriately 
generalized prescription. This can be accomplished by following through the arguments used 
in Q to derive (|50|): first of all the scattering considered below is diagonal in the sense that 
whenever a spinon or holon bounces off a wall, it merely changes its rapidity A to —A and 
acquires a phase-shift. Secondly the wave-functions of the excitations on our finite interval 
with fixed boundary conditions are standing waves of states with opposite rapidities, which 
leads to the condition 

exp(zLp(Ai))5 , (A a , A 2 )K 1 (A 1 , h x ) = exp(-tLp(X 1 ))S(~X ll X 2 )K 1 (-X 1 , h L ) , (62) 

where p(X) is again the physical momentum of the corresponding (infinite) periodic system 
and S(X, v) are the bulk scattering matrices. The extraction of the boundary phase-shifts from 
evaluating the counting functions at the spectral parameters of the holons/spinons is much 
more problematic than for the half-filled case as now the integration boundary B corresponding 
to the charge Fermi sea is finite and the issue of how to treat the corrections to B arises. 



V. CONCLUSION 

In this work we have studied zero-temperature boundary effects in an open supersymmetric 
t-J chain with boundary fields. Surface contributions to ground state properties were evalu- 
ated as functions of the boundary fields and the phase-shifts acquired by holons and spinons 
scattering off a boundary were determined. It also would be interesting to extend the analysis 
to finite temperatures. This appears to be difficult as among other things the usual expression 
for the entropy [^] has to be modified in order to project out the spurious states. The best 
way to deal with these problems seems to be a Thermal Bethe Ansatz analysis |I0[ 



The (subleading) finite-size corrections to the ground state energy, and energies of low- 
lying excited states, can be evaluated by means of the Euler-Maclaurin sum formula as well. 
This was done for the case of the aa-boundary conditions and zero bulk magnetic field in 
|4li . That analysis, which can be straightforwardly extended to the other sets of integrable 
boundary conditions considered here, allows to make contact with Conformal Field Theory in 
a geometry with boundary [42 . 
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APPENDIX A 



In this appendix we derive the expression (|T3p for the ground state energy. The ground 
state for given magnetic field H and chemical potential p is obtained by filling all vacancies 
for the integers I* from 1 to J max = N± + Nh and all vacancies for the integers J 7 between 1 
and J max = Nh. Inserting this description into the Bethe equations (0) and then subtracting 
susequent equations for a and a + 1 and 7 and 7 + 1 we obtain the following equations for 
the densities g s (X a ) = jjt — 1 and g e (X^p) 



L(\ a+1 -x a ) wy, 1 l ( a« 1 -aV ) ) 
g s (X a ) = 2a 1 (A Qt ) - -^a 2 (A Q - A^) + a 2 (A a + \p) + — ^a x (\ a - A^) + a x (A Q + A^) 




QcVW) = lE ^ - A «) + + A «) + ^Z^'^^) • ( Al ) 

Here k^-(A) and u^-(A) are the derivatives of >4j ( A ) and ^(A) defined in ([TT|) respectively. 
Now we follow [ffl and rewrite (IAU) in terms of a set of "doubled" variables 



-X Nl+Hh -a OL = 0, . . . , N t + N h 

Xa-^-N, a = N l + N h + l,..., 2(N l + N h ) 

-A$ ^ 7 = 0,...JV A 
v?> = { h " 7 (A2) 



A 7 -jv h 7 = N h + 1, . . . 2N h , 

where we defined Ao = and Aq 1 '* = 0. Now we take the thermodynamic limit of the equations 
( | Al| ) written in the new variables. This is done by using the Euler-Maclaurin sum formula 
to turn sums into integrals (see e.g [(43|,[y|(44]]). Care has to be exercised in order to take 
into account the fact that terms depending on the spectral parameters located at zero must 
be subtracted explicitly. After some manipulations we arrive at following coupled integral 
equations for the densities p s {y a ) = j- ( — ~ — \ an d p c (^ 1 ^ 



p s (X) = 2ai(A) - / dp a 2 (\ - p) p s (p) + / dp ai(A - p) p c (ji) + T I h a 2 (A) 

-A+ J-B+ L \ Z71 



Pc(X) 



f\p ttl (A - p) p s (p) + \ - 01(A) J , (A3) 

where v4 + and B + are the spectral parameters corresponding to the maximal taken integers 
I a and J 7 plus |. Higher order terms in the Euler-Maclaurin expansion have been dropped as 
they turn out to not contribute to the surface energy. 

The energy Der site (|6[) in the thermodynamic limit can be expressed in terms of the 
densities p SjC asM 



**We again turn the sums into integrals by means of the Euler-Maclaurin formula. 
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l = ~ J-A+ ^ ai ^ v > ~ ps ^ v > + J_ B+ dv ~ ip 

1 ' + + 2 -//--. (A4) 



LI 2 4 

Note that we have divided the bare energies by two due to the fact that we are working with 
the densities of the doubled variables. We also explicitly subtracted a term + ^ — 2) to 
compensate for the spurious roots at zero. To proceed further we rewrite (|A3|) as an operator 
equation 



~K SS -K sc \JPs\_ M° } 

-%s i-kJ*\ p J U 0) 



(A5) 



where I is the identity operator and * denotes convolution with the appropriate kernel e.g. 
-K ss * p s 



f A ^+ dp a 2 (A — p) p s (p) and where 



,<»> = 2 Bl (A) + i + a 2 (A)) , ,f = I - a l( A)) . (A6) 

We note that the above integral kernels are all symmetric. In what follows we use the 
shorthand notation (id — K) a bpb = p^ for (|A5|) . Let us now define quantities e s (X) and 
£ C (A) through the integral equations (id — K) a bEb = 4 0) , where e^(A) = -27rai(A) + H and 
£( c\^) = A* — 2" • With x — ~2 ~ T + -^v + ^ the energy per site can now be written as 



T = | 51 / \ 4 O) (a0p&(a0 - A* - y + | 

b=s,c ~~ 

= \ E !°\dpet\p)(id-K)^^p^ 

a,b=s,c b 



H y 



~ £ ^ ^)pi 0) (A*) - A* - f + 7 , (A7) 



where = A + and = B + . In the thermodynamic limit the ground state is obtained 
by minimizing E with respect to the integration boundaries A + and B + (see e.g. f44|), i.e. 

-£^p = 0, where C s = A, C c = B. From this fact it follows that the integration 

boundaries in ( |A7| ) can be replaced by C a with error of 0(L~ 2 ), which does not affect the 
surface energy we are after. In other words we can replace by C a (up to 0(L~ 2 )) due to 
the fact that the dressed energies vanish at the integration boundaries. This finally establishes 



APPENDIX B 

In this appendix we outline how to evaluate the integral fl38|). We first note that via Fourier 
transform the following equality can be established 
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dfi a s (n) e 8 (n) = -2nG s {0) + — + 2 / dp [G 1+S (fi + A) - a s (fi + A)] y(fi) 
-A 2 Jo 

+ [ B G s (jm) e e (jM) . (Bl) 

J-B 

The last term on the r.h.s. is readily evaluated by using the fact that Gs{fJ>) is smooth around 
zero and thus can be Taylor expanded 



/ 



B G s (fx) e c (») = -2C(3)5 3 G 5 (0) . (B2) 

B 



The second term on the r.h.s. (which will be denoted by R 2 in the following) is more difficult 
to treat. In the region where S 1 < i we can use the asymptotic expansions G 1+ s{n + A) ~ 
47r(^+% 2 an d a s(r l + A) ~ 2-K(f+A) 2 ^° determine the leading contribution to the integral. We 
then Laplace transform jjj^f an< ^ & fter some manipulations arrive at 

roo I roo 

J (A+fj) 2 y ^ = Jo dx xe • ( B3 ) 

Due to the strongly decaying factor e~ Ax the leading contribution to this integral clearly comes 
from the region around x — 0. Expanding y + (ix) around x = we finally obtain 

i?2 = — - — ; — — + subleading terms . (B4) 
2 In H 

In the region S 3> A the above strategy for determining R 2 is not applicable. Instead we 
Fourier transform and arrive at 

(B5) 



We again have a strongly decaying factor in the integrand (recall that 4^- ^> 1) so that we can 
expand the other terms around u = 0. The problem we run into now is that the subleading 
terms in the expansion of y(A) contribute in an important way. The leading contribution 
of yi to i?2 is found to be (after expanding yf around cu = only elementary integrals are 
encountered) 

H Hln(S-l) 2H\nH 

The second term is precisely cancelled by the leading contribtuion from y 2 to R 2 , whereas the 
further terms are all small compared to — |Jpzf) • 
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